Gluon saturation and inclusive production at low transverse momenta 
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In this letter we suggest the generalization of fcT-factorization formula for inclusive gluon produc- 
tion for the dense-dense parton system scattering. It turnes out that the soft gluon production with 
transverse momentum pr is suppressed by additional Sudakov-like factor that depends on Px/Qs 
ratio in a good agreement with the first numerical calculation in Colour Glass Condensate approach 



, by J. P. Blaizot, T. Lappi and Y. Mehtar-Tan. 



It is well known that our approach to inclusive production of an gluon jet is based on kx factorization[l|-[4| which 
leads to 

da 27ra s 



^ dyd2pT p 2 ± j Jk T 4 1 kr) ^ (x2-,pt - hr) (1) 

Q-T h • . 

where <Pq are the probability to find a gluon that carries Xi fraction of energy with k± transverse momentum and 
a s = a s N c /ir with the number of colours equals N c . 

In the framework of high density QCDjH- 13 1 the kx -factorization has been proven [l4| (see also Refs. 15, 17-2(|) 



for the scattering of the diluted system of partons, say for virtual photon, with the dense one. Such scattering is 
characterized by two scale of hardness: the saturation momentum of the dense system Q s and the px of the produced 
gluon. The dense-dense parton system scattering has three scales of hardness: two saturation momenta and px] and 
the kr -factorization has not been proven for this process. The most dangerous region is for px smaller than both 
£C) \ saturation momenta ( p t < Q\, s < Q 2 s) where we did not expect that kx factorization will work. However, for 
Qi.s < Pt < Q2s we are dealing with scattering with two scales of hardness and we can expect that kx factorization 
is valid here. In this paper we address this problem and suggests the generalization of the kx factorization (see below 
Eq. Q) for p t < Qi, s < Q 2s - 

As it was shown in Ref.fbl] <f>Q (xi;k±\ can be written through dipole scattering amplitude N (xi,rj_:b), where 

r±_ is the dipole size and b is the impact parameter of the scattering. This relation reads as follows 



X 

where 



0&(xi;fcr) = J d 2 bd 2 r ± e^ T ^ W 2 ± N^ (y t ^ln(l/x t );r ± ;b) (2) 



N%( yi = ln(l/xi);r ± ;b) = 2N ( yi = ln(l/x i );r x ;b) - N 2 (y 4 = \n(l/ Xl ); r±; b) (3) 



N (yi = \n(l/xi); r±; b) is the dipole -hadron (hi) scattering amplitude which satisfy the Balitsky-Kovchegov equation. 
Using that N (xi ,r±;b) is a function of r\ we can rewrite Eq. ([2]) in the form 



iU;k T ) = ^— [ d 2 bd 2 r ± e^-^—J-r x J-N^(y i = ln(l/x i );r ± ;b); (4) 
V / a s 47T J r± or± or± 



irj_ — oo 



d 



- ik T I dr x d<l>cos ( j>e ik ^ ms *r ± —N%(y i =\n(l/x i y,r ± -,b)y, (5) 

At r± — > the amplitude Nq approaches the solution of the DGLAP equation which in our case corresponds the 
double log limit of the BFKL equation. Therefore 

d 2 br ± J-N%(y i =]n(l/x i );r±;b) -> DLA (6) 
or _l 



cc r ± -^- exp( v /4a s ln(l/ a ;) / {r\\ 2 QC D ) - ln(l/(ri A£ CZJ )) ^> 
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FIG. 1: Function cj>o versus transverse momentum in 
different saturation models: GBW denotes the Golec- 
Biernat and Wusthoff model (see Ref.JIl]) while CGC de- 
notes the model suggested in Ref . [23l] . These two model 
have different values of the saturation momentum and 
the picture illustrates that 4>{x,pt) has maximum at 

PT = Qs- 



It is worth to mention that the impact parameter dependence enters in Eq. ([6]) as a factor and cannot change our 
claim that this term vanishes at rj_ — > 0. 

At large r± 1 - N G cx exp f - Cln 2 (r 2 Q 2 ) ) (see Ref.[IH) and, therefore, 



d 

Finally we see that the first term in Eq. ([5]) vanishes and we have 

d 



d2br ±^~ N G (tt = In(l/Jij)in;i) (7) 



<f>% (xn k T ) = k T dr ± Ji {k T r±) r±—N& { Vi = \n{l/x l );r ± ;b) (8) 



dr 

From Eq. ([8]) one can see that at kx — >• <po <x k\ while at large values of kr 4>g °c l/k T — > 0. Such behaviour 
of <Pg means that it has maximum at kx ~ Q s - The numerical calculations [33j| (see Fig. [T]) confirm this claim. 

Having this feature of <pG in mind we see that at pr <C Q s Eq. (fTJ) gives 

da 2ira 



dy d?p T p\ 



d 2 k T 



^ (x.-kr) ^ (x 2 - k T ) (9) 



One can see that at px — > the cross section tends to infinity. Since we are talking about inclusive cross section 
generally speaking such situation is possible and it corresponds to increasing multiplicity of soft gluons. However, in 
framework of gluon saturation it looks strange. As we have discussed above, the main contribution to <j>c give the 
gluons with transverse momenta of about Q s while the gluons with small values of kx are suppressed. In other words, 
the correlation length between emitted gluon is of the order of l/Q s and we expect that emission of gluons with the 
wave length larger that l/Q s should be suppressed. Of course, soft gluons with px <C Q s could be emitted in the 
final state but they will not propagate through the medium since the cross section is large (aT (b) cx Q 2 s /p T 3> 1). 

In this letter we will show that simple formula of Eq. ((T|) should be changed and a new double log suppression factor 
(T) should be added. Therefore, the inclusive cross section has a form (see Fig. [2]-a) 

= -^r 1 J d2( li d 2 ?2<K<?i -<h ~Pt) 4>g {xx;qx) 4>q (x 2 ; q 2 ) T (^) ( 10 ) 



dy d 2 p> 



In Eq. (|10p we assume that Q s ,i ~ Qs.2- The appearance of T in inclusive production was found in 1980's [2J426J 
and it is related to the fact that the emission of some gluons is suppressed in the process. In our case the emission of 
gluons is suppressed with the value of the transverse momenta ( pr,i) in the region: px < pr,i < Qs (see Fig. [2}b, 
where the gluons which emission is suppressed, are denoted by the dashed lines). Actually, the emission of gluons 
with small values of px,i has been taken into account in functions <pQ but they result in suppression of the emission 
for such gluons and we do not need to account separately for them. 
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FIG. 2: Inclusive cross section: <fi — <j2 = Pt- 



Eq. (jlOp says that the emission of the gluon with px < Q s is suppressed and only gluons with p t > Q s gives the 
contribution to the inclusive production. For such gluons the kr factorization works. This qualitative features of 
Eq. (I10[) has been confirmed by the first numerical calculculation that found the deviation from kr factorizxation (see 
Ref.[28|). These calculation shows that for pr < Qs the gluon production is suppressed while for px > Qs the kr 
factorization works perfectly well. 

We calculate the first diagrams of Fig. [2}c to illustrate the double log contribution. 

This diagram is equal to 



A(Fig.m-c) = -T^r-AUv^U^afa 

(Z7T) 4 I V 



Nr. 



fa- 



be 



dk+ dk~ d 2 k 7 
2 '""V J Iqt -<h 
^ r 7: q lP q 2p 



(11) 



(k+k- - fcf, - ie) - fe+)( ? r - k-) - (ft - k)\ - ie) - fc+)(<fc - k~) - (q 2 - kf T - ie) 



In Eq. (Ilip we used that at high energies the propagators of gluons with momenta q\ and qi can be written in the 
form % [23 



(qi) = - 4. 72—; (92) - 



q-2 qt 



q 2 qj q ltT 

In leading log(l/x) approximation we have the following kinematic constraints: 

qt > qt; ql > it <h = pt; qiql < ?i,t; ^ ?2 < <?2,t; 



(12) 



(13) 



Integration for kr > Q s leads to renormalization of the coupling QCD constant. Therefore, we are interested in 
the kinematic region where 

q it T « Qs > k T > p T (14) 
Having in mind Eq. (fT"3")) and Eq. (fT4")l we can rewrite Eq. (fTTj) using 

( qi - kf +ie= (k- - q-) (k+ - q+) - Q 2 S - ^e; (q 2 - kf + ie = (k+ - q^)(k~ - q 2 ) - Q 2 S - ie; (15) 

One can see that for < k~ < q 2 the poles in fc + in Eq. (jlip are situated in diffrcnt semiplanes and we can 
close the contour in complex plane k + on the pole from (92 — k) 2 = 0. It gives 

QI QI 



k + 



k - q 2 



(16) 



g 3 n N c 



(2tt) s 



fa 



be 



dkrp elk 



2 Prl2 ( q +k- + QI) (kr&. + k, 



(17) 
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One can see that from the kinematic region given by 

kj, q 2 



n , » k~ » Qi/gf; Q 2 » k 2 T » p£ (18) 
we have a double log contribution, namely, 

A (Fig. m-c) = - ^3 ^ / a6c 4p2 qf^a.p I>, 7il/ r 7 , a , p ^ In 2 (Q 2 /p|) (19) 

Direct calculation of sum over gluon polarization in Eq. (|T51) leads to 



^-^ r p,a,/3 r /3,7^ r 7,a,p92,p = 4 ^ 2 ) l\, (20) 



where 

r* = 2 fe.i - J^) (21) 

is famous BFKL vertex [27} . 

Collecting all factors and using the notation T° bc = 2g f a b c F v we obtain for the diagram of Fig. [2}c the following 
expression 

A (Fig. m-c) = Tf c In 2 (q 2 /^) (22) 



where a s — a s N c /n. Using the well known technique (see Refs.[2J, |25|) we obtain 

T (|) - M-y^ («/«*)} < 23 > 

This result follows directly from the generalization of Low theorem for soft photon [2{| for high energy scattering 
(Gribov's theorem 30]). It says that if px of emitted photon smaller than any typical transverse moneta in the process 
the cross section of emitted photon is equal to 

a du cPpr 

ct 7 = — a (24) 

2ir u> p T 

where erg is the cross section for the process without photon. This theorem has been generalized to emission of 



gluons (see Ref. |3l|). In our case the typical momentum scales of the process are Q s ,i and Q s ,2 and, therefore, 
gluons with p? <SC Q s ,i and Q s 2 are emitted independently according to the Poisson distribution. The emission of 

one gluon will be suppressed by exp(— (n)) where (n) = ^ In 2 [Q^/p T j is the average number of emitted gluons. 

In the case of two different saturations scales Q s ,i and Q s .2 the average multiplicity (n) = hi 2 ^Q 2 m j n /PxJ 

where Q s ,min = niin{Q s .i, Q s .2\- For dilute-dense scattering the value of Q s ,min vanishes and we have no additional 
suppression. This feature is related to the fact that the BFKL evolution has two branches one of which leads to 
decrease of the typical transverse momenta of gluons. 

In conclusions, we would like to summarize that we suggest Eq. (|10p which violates the kx factorization for kx < 
Q s ,i andQ S; 2 but it is in a perfect agreement with the numerical solution for the inclusive production in Colour Glass 
Condensate (28|( see Fig. [3]). We would like to emphasize that our result is based on general grounds in QCD and 
reflects the fact that emission of soft gluons with transverse momenta Q s > pi y T > Pt has been taken into account 
in functions 4>g (see Eq. (|10I0 and it should not be included again in fc^-factorization formula. 
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FIG. 3: The single inclusive cross section 
p 2 T dN/d 2 Pt versus pr- The low curve shows the 
exact solution given in Ref. [2^] . The upper 
curve presents the result for the inclusive cross 
sections from fcT-factorization. Both these curves 
are taken from Ref. [28I ] and we are very thank- 
ful to Tuomas Lappi who shares with us the data 
for these curves. The blue curve (the shortest 
one) shows the result for Eq. f| 10|) assuming that 
a s = 0.1 and Q s — 1.5g 2 fi. 
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